Entanglement of nanomechanical oscillators and two-mode fields induced by atomic 

coherence 
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We propose a scheme via three-level cascade atoms to entangle two optomechanical oscillators as 
well as two-mode fields. We show that two movable mirrors and two-mode fields can be entangled 
even for bad cavity limit. We also study entanglement of the output two-mode fields in frequency 
domain. The results show that the frequency of the mirror oscillation and the injected atomic 
coherence affect the output entanglement of the two-mode fields. 
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I. INTRODUCTION 



Entanglement among mesoscopic even macroscopic 
systems have always been an attractive topic since the 
birth of quantum mechanics. Schrodinger's cat para- 
dox is a well known example of it (entanglement between 
macroscopic cat state and microscopic state) . Cavity op- 
tomechanical system is an important candidate for the 
study of quantum mechanical features at mesoscopic even 
macroscopic scales. Recently, there has been considerable 
interest in studying entanglement in mesoscopic systems 
[1-9]. Refs. [3, |4| investigated the optomechanical en- 
tanglement between the field of an optical cavity and 
a vibrating cavity end-mirror. Many proposes were put 
forward to entangling two mirrors in a ring cavity [6| , en- 
tangling two mirrors of two independent optical cavities 
driven by a pair of entangled light beams [7[ , entangling 
two mirrors of a linear cavity driven by a classical laser 
field Q, entangling two separated nanomechanical oscil- 
lators by injecting broad band squeezed vacuum light and 
laser light into the ring cavity Q ■ Most recently, Ref . [|[ 
observed strong coupling between a micromechanical res- 
onator and an optical cavity field, which could push for- 
ward the possibility of entangling two mirrors via optical 
cavity field. 

On the other hand, atomic medium play important role 
in the interaction of cavity electrodynamics. When the 
atomic medium is trapped in a cavity with micromechan- 
ical oscillator, some interesting phe nomena have been 
revealed |10l - [l3j . Genes et al [10] suggested cooling 
ground-state of micromechanical oscillators by resonant 
coupling of the mirror vibrations to a two-level atomic 
bath. Ian [ll[ found that two-level atoms effectively 
enhance the radiation pressure of the cavity field upon 
the oscillating mirror. Most recently, Hammerer et al 
[l2| and Wallquist [HI have shown the possibility to 
achieve strong cavity-mediated coupling between a sin- 
gle microscopic atom and a macroscopic mechanical os- 
cillator. In addition, Genes 14] studied tripartite entan- 
glement among atom, mirror and cavity field. As we all 
know, atomic coherence results in many interesting phe- 



nomena |l5[ such as electromagnetically induced trans- 
parency (EIT), correlation emission laser (CEL), laser 
without inversion (LWI). A lot of works have been done 
on the entanglement directly induced by atomic coher- 
ence (l6l - [20| . When one of the mirrors of the cavity is 
movable, the atomic coherence effects are not studied be- 
fore. In this paper, we propose a method to entangle 
two macroscopic mirrors via microscopic atomic coher- 
ence. When atomic beam with cascade configuration is 
injected into the two-mode cavity, the two-mode fields as 
well as th e op tomechanical oscillator are entangled. In 
Refs. [Ill, Il3| . authors treat the cavity field as a quan- 
tum bus and give an effective coupling between oscillator 
and the single atomic motion. Instead, we directly treat 
hybrid system. Our study show that the initial atomic 
coherence and the frequency of the mirror's motion affect 
the entanglement of the output fields. 



II. MODEL AND THE HAMILTONIAN OF THE 
SYSTEM 



The system under study is a two mode cavity with one 
fixed partially transmitting mirror and two movable per- 
fectly reflecting mirrors, sketched in Fig. 1. The atomic 
medium with cascade configuration is injected into the 
cavity and interacts with two-mode cavity fields with de- 
tuning Aj, respectively (i — 1,2). The Hamiltonian of 
the hybrid system reads 



H = 



J=l,2 .7=1,2 



(i) 



J=l,2 J=l,2 

E hEiUu + h(gia ba a\ + g2^cb<A + h.c.) 



i—a,b,c 



The first term describes the energy of the two cav- 
ity modes, with lowering operator aj, cavity frequency 
LJj, and the decay rate Kj. The second term repre- 
sents the energy of the two mechanical oscillators at 
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frequency u> mj , and Pj and Qj are their position and 
momentum operators. The third term describe the two 
driving laser with frequency u>l 1 and wi 2 , respectively. 
The forth is the radiation-pressure coupling with rate 
Xj = T7^/mJ~~"' an< ^ l as * summation describes the en- 
ergy of the atoms, and the last term is interaction be- 
tween the atom and the cavity fields, where au = 
is the spin operator of the atom. 

In interaction picture, we have the Hamiltonian 



Hi = 



Mi 

+H(giab a a\ + g2U c ba\ + h.c.) 
+h(Si — Ai)a\a\ + fi(<5 2 — A 2 )a 2 a 2 
+ie\h(a\ — Oi) + i£ 2 7j(a 2 — a 2 ) 



(2) 



(Pi 



2 v 1 ^' 2 
+hXiQia\a 1 + hx2Q2a\a 2 ], 

where A x = E a - E b - a>i, A 2 = E b - E c - cj 2 , <5i = 
E a — Eb — u>l 1 , 8 2 = Eb — E c — ujl 2 



Thus Sj — Aj = 



ujj—ujLj means the detuning between the classical driving 
fields and the cavity field. The dynamics of the system is 
determined by the following quantum Langevin equations 



Qi 
Q2 
Pi 
Pi 
o'i 

a" 2 
0"cb 



,Pl 



(3) 



6, 



?1 - Tm-Pl + 

-\ki + ?Ai]ai - igiaba + £\ 
-[k 2 + ?A 2 ]a 2 - i.g 2 cr cfc + e 2 
-(7 + iSi)a ba - ig\ai{a b b - cr aa ) + J.g 2 a 2 cr ca , 
-(7 + i<5 2 )er c 6 - igia\a ca - 15202 (o- cc - &bb), 



'2K 2 a 2m , 
,t. 



where Aj = Sj — Aj + XjQj (J — 1, 2). The quantum 
Brownian noise £1 and £ 2 are from the coupling of the 
movable mirrors to their own environment. They are 
mutually independent with zero mean values and have 
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FIG. 1: The sketch of the system and the atomic configura- 
tion. 



the following correlation function at temperature T 

Sjklm 



&(*)&(*')> 



0J„ 

j,k = 1,2. 



Z7T AKbJ- 



(4) 



The two cavity modes decay at the rate K\ and k 2 , and 
aim (a2m) is the vacuum radiation input noise, whose 
correlation functions are given by 



(<xt. n (t)a jin (t')) = NS(t-t'), 
(a jm (t)aUt')) = (N+l)S(t-t% 



(5) 



where N = {exp{huj c /k B T) - l] -1 . 

In order to obtain the steady state solution of Eq. ((3|) , 
we calculate the last two equations to the first order in 
Qi (i — 1,2), i.e., using linear approximation theory [ToT ]. 
which means that in the last equation of for the 
terms that the a\j multiply a (aJ), we use zero order 
(<Jij) substitute of ay. We assume that the atoms are 
injected into the cavity with the state p a — p° aa \a){a\ + 
Pcc\ c )( c \ + Pca(\ c )( a \ + l a )( c l) at injection rate r a . The 
last two equation of ([3]) can be rewritten as 

o"6a = -(7 + *<5i)cr 6a + ig\r a p\ a ai + ig 2 r a p ca al, (6) 



&cb 



-(7 + *<5 2 )cr c6 - igir a p Q ca a\ - ig 2 raP u cc a 2 



By combining it with Eq. (|3|) , we finally have the steady- 
state mean values of the system as 



PI = 0,P| = 



, Q2 



-X2\a s 2 ? 



(J) 



s 2c ei + e 2 ui 

U1U2 + SlaS 2c ' 

s* la e 2 - e\u2 

U\U2 + S2cS* a 



with 



u, 



Sla 



• s 2r 



(0) 

gig2r a pca 
7 + iSi 

Ki + iAi - 
k 2 + iA 2 H 



,1 = 1,2 

a 2 r o (0) 
7 + i5i 

2 (0) 
ff 2 r aPcc. 

7 + i&2 



(8) 



From the expressions of (7) and (8), we see that if g\ — 
g 2 = (no atoms within the cavity), the steady value 
of a,j (j = 1,2) will have the same form with [3, El, he., 

If the injected atoms have no coherence 



3 k j -H Aj- 



between their levels, that is to say, p^ia = (ui = 0). 
By writing the last terms of s la and s 2c into real and 
imaginary parts, we know that the existence of atomic 
medium only affect the effective decay rate of the photon 
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and the radiation pressure of the cavity field upon the two 
mirrors. However, if the injected atoms have coherence 
between the level \a) and |c), i.e., p^a 7^ 0, the two mode 
fields will be correlated so that the two mirrors will be 
dependent each other, see (7) and (8). 

We write each operators of the system as the sum 
of their steady-state mean value and a small fluctu- 
ation with zero mean value. The fluctuations can 
be calculated analytically by usin g t he linearization 
approach of quantum optics [H, HH, [2l[, provided 
that the nonlinear effect between the cavity field and 



the movable mirrors is weak. We can linearize the 
fluctuations around the steady state using Eqs.Q 
and obtain a set of linear quantum Langevin equa- 
tions for the fluctuation operators. We define / = 
{5Qi,5Pi8Q 2 , SP 2 , S Xi ,6Y 1 ,6X 2 ,SY 2 ,SU 1 ,6U 2 ,SV 1 ,6V 2 ) T 
A set of linear quantum Langevin can be written as 

f = Af + B (9) 

where 



A = 
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B = (0, £i,0, £2, \f2k~l5X lin: s /2k~l8Y lin , y/^SX^, 
^SY 2in , 0, 0, 0, 0) T . We have defined 



Entanglement of the nanomechanical oscillators 
and the two-mode fields within the cavity 



X 3 = -L( 0i + oj), Y j = -^ ] (a j -a]),j = 1,2,(10) 



V2i 



Vi 



1 1 

—r^K^ba + (Tab), U 2 — —=-{(Tba ~ & ab) , 

V 2 V 2i 

1 1 

—7^{ (7 ba + &ab),V 2 — —=-{(T c b — (Jbc)- 

V 2 v 2i 



A is the drift matrix. The system is stable only if the real 
part of all the eigenvalues of the matrix A are negative, 
which is also the requirement of the validity of the lin- 
earization method. Because of the atomic cascade form, 
the system can be a amplifying system, that is to say, for 
some parameters, the eigenvalues of drift matrix A can 
have a positive real part. So, we must carefully eliminate 
the parameters region in order to retain it within the sta- 
bility conditions. Due to the twelve dimensions of drift 
matrix A, it is not easy to obtain the analytic expression 
of the requirement as Q eigenvalues. We will guarantee 
the requirement via numerical method. All the parame- 
ters chosen in this paper have been verified to satisfy the 
stability conditions. 



We investigate the nature of linear quantum correla- 
tions among fields and mirrors by considering the steady 
state of the correlation matrix of quantum fluctuations in 
this multipartite system. The quantum noises £ and a\ n 
are zero-mean quantum Gaussian noises and the dynam- 
ics has been linearized, as a consequence, the steady state 
of the system is a zero-mean multi-partite Gaussian state. 
We defined V^-(oo) = !K/i(oo)/j(oo) + /j(oo)/;(oo))], 
the element of covariance matrix. One can obtain vari- 
ous correlation information from it, and we will show the 
entanglement of the two mirrors as well as the two-mode 
fields after obtaining the covariance matrix. 

We now recall entanglement criteria of continuous vari- 
able proposed by Duan [22| and Simon (23[. According 
to [22] , a state is entangled if the summation of the quan- 
tum fluctuations in the two EPR-like operators X and Y 
satisfy the following inequality 



(AX) 2 + (AY) 2 < 2. 



(11) 



For the mechanical oscillator we define X m = Qi + 
Q 2 ,Y m = Pi — P 2 while for the two mode fields Xf — 
X x - X 2 ,Y f =Y t + Y 2 . The criterion in Eq.(ll) can be 
directly detected in experiment via homodyne measure- 
ments 24]. Simon's criterion has more direct relation 
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with covariance matrix proposed in [23| and developed 
in [25]. For a physical state, the covariance matrix must 
obey the Robertson-Schrodinger uncertainty principle, 



V+ -/3>0, 



(12) 



where /3 



when we de- 



JO 1 

o j\ Wlth J = [-1 

fine vector / = (Sqi, Spi, 5q2, 5p2) T for two-mode fields 
23]. If a state is separable, partial transpose matrix V 
(obtained from V just by taking pj in — pj ) still obey 
the inequality in (|12p . The inequality requires that all 
the symplectic eigenvalues of the transposed matrix are 
larger than 1/2 in terms of the definition of Eq.(10). The 
symplectic eigenvalues can be calculated from the square 

roots of the ordinary eigenvalues of — 

(pv) MM>- so, 

if the smallest eigenvalue is smaller than 1/2, the trans- 
posed mode is then inseparable. For the two-mode Gaus- 
sian states, the violation of the inequality is a sufficient 
and necessary condition for the existence of entanglement 
between the transposed mode and the remaining modes. 

Neglecting the frequency dependence as it was pointed 
in [J , under the condition of Markovian approximation, 
the frequency domain treatment i s eq uivalent to the time 
domain derivation considered in |l0j. When the stabil- 
ity conditions (all real parts of eigenvalues of matrix A 
are negative) are satisfied, the steady state correlation 
matrix of the quantum fluctuation meet the following 
Lyapunov equation [J] 



AV + VA T = -D, 



(13) 



where D = Diag[0,'j mi (2N + l),0,-f mi (2N + 1), k 1 (2N + 
l),Ki(2JV + l),K2(2iV+l), k 2 {2N+ 1), 0,0, 0,0]. 

For simplicity, we choose all parameters of the two 
mirrors and the two mode cavities are the same, for ex- 
ample 3i = g 2 , w mi =_w m2 , uj Li = w l2J etc.. Then 
the effective detuning Ai = A2 = A. In Fig. 2, we 
choose the atoms injected into the cavity with initial state 
|^a(0)) = "^f (l a ) + l c )) which is the best coherence of the 
two levels. We numerically show the entanglement of the 
movable mirrors as well as the two mode fields induced by 
atomic coherence. Fig. 2a show that with the increasing 
of the coupling between cavity and the atoms, the deep 
squeezing in fluctuation move towards the side of small 
value of A. As already shown in 0, HI], many features 
of the radiation-pressure interaction in the cavity can be 
understood by considering that the driving laser light is 
scattered by the vibrating cavity boundary mostly at the 
first Stokes ujq — u) m and anti-Stokes loq +uj m . Therefore, 
the optomechanical interaction will be enhanced when 
the cavity is resonant with one of the two sidebands, i.e., 
when A = uj m , where entanglement between cavity field 
and mirror is enhanced. Here, due to the atomic medium 
existence, the two sideband move towards its middle. In 
addition, within the region of parameter in Fig. 2a we can 
observe that the larger value of the coupling, the deeper 



squeezing in the fluctuation. Undoubtedly, the entangle- 
ment is resulted from the interaction between atoms and 
the fields as well as the atomic level coherence, therefore, 
the larger coupling gi the deeper squeezing within our 
parameters. The enhancement of gi can be substituted 
by increasing the injection rate of the atoms. Fig. 2b 
show that the two-mode fields are entangled for nearly 
the whole region of detuning A (from to 1). Using the 
same group of parameters, from the solution of Eq.(13) 
we reconstruct two-partite covariance matrixes for mir- 
rors and two-mode fields and show their smallest eigen- 
values under partial transpose, shown in Fig. 3. One can 
see that Fig. 3 are almost the same in shape with Fig. 2, 
that is to say, the entanglement judged with two criteria 
are equal. By showing the entanglement with two crite- 
ria, we prove our correction of calculation and convince 
the entanglement existence. 

For cascade atomic configuration system, the pro- 
cess of radiation two-mode fields are correspondence to 
parametric-down-conversion process. So, one can easy 
understand the two- mode entanglement. For fixed mir- 
ror, it has been known that two cavity modes are entan- 
gled when they interact with cascade three-level atoms 
[16j . Here, we show further that when one mirror is mov- 
able, the two-mode entangled continuous-variable state 
still hold. In our numerical simulation, we find that there 
is a threshold coupling between atoms and fields (un- 
der these group of parameters), below that we have no 
two entangled mirrors, which means that the appropri- 
ate coupling between atoms and cavity is favorable if we 
would have entangled mirrors. Thus, the radiation fields 
emitted from cascade atoms can be entangled, and the 




0.6 0.8 
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FIG. 2: (Color online) (a) shows the entanglement of the 
mechanical mirror where gi = 2ir x 1.5 x 10 5 Hz (solid), 2ir x 
1.7 x W 5 Hz (dash-dot), 2tt x 2.0 x 10 5 Hz (dashed), (b) plot 
the entanglement of the two mode fields where gl — 2tt x 2.2 x 
10 5 Hz . For all the plots, the other parameters are L = 5mm, 
m = 20ng, m = K2 = 2tt x 215KHz, tj mi — cj m2 = 2ir x 
10MHz. The wavelength of the laser A = 810nm with power 
lOmW, and the mechanical quality factor Ql — — 6700, 
ra = 2000, r = 12MHz, 81 = S 2 = 4MHz, T = A^JiK . 
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entanglement of the two-mode field can be transferred 
into the mechanical mirror because of the radiation- 
pressure interaction. As a consequence, we have both 
the entanglement of two-mode fields and two mirrors. 



B. The output two- mode field entanglement with 
optomechanical oscillator 

Because one can not direct access to the intracavity 
fields, if one would like to use the entanglement, he need 
to export the entangled fiedls first; thus only the output 
entanglement of the optical cavity has practical meaning 
. In addition, by means of spectral filters, one can always 
select many different traveling output modes originating 
from a single intracavity mode, which offer us the oppor- 
tunity to easily produce and manipulate a multipartite 
system, eventually possessing multipartite entanglement 
4]. So, in this section we study the entanglement of the 
output field. The input-output relation for the two mode 
fields are 

Sajout = \/2KjSdj - 5a Jln ,j = 1,2. (14) 

The linear equation of © can be solved in the fre- 
quency domain by Fourier transformation with the so- 
lution f{oS) ~ {—iuj — A)~ 1 B. In the interaction picture 
co represents the detuning from the cavity frequency [29| . 
Considering the relation of (|14|) . we have 

f out (uj) = C{-iuj - AY 1 B - E, (15) 

where C = diag(0, 0, 0, 0, v 7 ^, V^i, V^2, V^2, 0, 0, 0, 0) 
E = [0, 0, 0, 0, SX lin (oj),SY Un {oj),dX 2in (uj), SY 2m (uj), 0, 0, 
0, 0] T . We can finally obtain the output correlation ma- 

trix v° ut (uj) = iK/r'H/^v) + /rV)/r*M>], 



0.5 




FIG. 3: (Color online) The entanglement of the mechanical 
mirrors and two-mode fields via the smallest eigenvalues of 
partial transposed matrix. All the lines and parameters are 
correspond to Fig. 2. 



and the squeezing spectrum 

S OU t(u) = ^[5X f (uj)6X f (oj') + SX f (oj')5X f {uj) 

+5Y f (oj)5Y f (oj') + SY f (u')5Y f (cj')} (16) 

can be calculated from the correlation matrix. The spec- 
trum is defined in a frame of 

Fig. 4 shows the two-mode output squeezing spectrum, 
which can also be considered as a measure of entangle- 
ment in frequency domain [l^]. In Fig. 4a, we notice that 
deep squeezing (usually means that the maximum en- 
tanglement) happen when w = 0, that is to say, when 
the output frequency is equal to that of cavity frequency, 
we have maximum entanglement. It is very interesting 
to observe that the deep squeezing are affected by the 
frequency of optomechanical oscillator. The larger fre- 
quency of the mirror, the smaller entanglement. We can 
deduce that if the mirrors is fixed, we should have the 
deepest squeezing in fluctuation of the output fields. The 
more higher frequency of the mirrors, the more entangle- 
ment can be transferred from the fields to the mirrors 
so that the smaller entanglement is leaved in the output 
fields. 

Fig. 4b shows the squeezing spectrum affected by the 
initial atomic state. As we have discuss above, the entan- 
glement of two-mode fields as well as the two mirrors orig- 
inate from the atomic coherence. Undoubtedly, the ini- 
tial atomic coherence affects both the intracavity entan- 
glement and the output entanglement. For the best co- 
herence of atomic initial state p° a = ^(\a) + \c))((a\ + (c|), 
we have one valley of deep squeezing when the out- 
put frequency is the same as to that of cavity (w = 




1 -6\i o o.i 14 -o.i o o.i 

co/co co/co 



FIG. 4: (Color online) (a) polts the output squeezing spec- 
trum for several values of frequency of the mirrors where 
w mi = = 2tt x WMHZ (solid) 2iv x WMHZ (dashed), 
(b) shows the squeezing spectrum under different input- 
atomic state . For solid line p aa = Pec = pac = 0.5, 
for dashed line p aa = |, p cc = |, p ac = |. In Fig.4b, 
cj mi = uj m2 — 2tt x WMHZ. For all the line in Fig. 4, 
gl = g 2 = 2n x 2.0 x 10 5 Hz , A = 0.8cj mi , and the other 
parameters are the same with Fig. 2. 
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in Fig. 4a). But for the input atomic state p° a = 
i|a)(a| + |(|a)(c| + |c)(a|) + ||c)(c|, the maximum squeez- 
ing have two valley. As to the split of the squeezing spec- 
trum, we can understand it from initial atomic state. The 
initial population in \a) and |c) are no longer equal so 
that the average of steady state \Xf\ (|l^ s |) for the two 
modes are not yet equal. Then, the fluctuation the 8X\ 
(SYi) and 6X2 (SY2) can be different so that the squeezing 
spectrum of the product of the fluctuation change from 
one degeneracy valley into two maximum squeezing. The 
split of the squeezing spectrum from one valley into two 
minima also have been shown in the case of asymmet- 
ric loss for each mode, i.e., K\ ^ K2 Our group 
had revealed the similar split resulting from asymmetric 
detuning and asymmetric atomic initial state [19| . 

During above numerical simulation, we choose param- 
eters of the system based on recent experiment for op- 
tomechanical system [j| , also considering the parameters 
used in theoretical papers [J, |9| . As to the coupling be- 
tween the atoms and the cavity, we use g = gl = gi = 
2tt x 2.8 x 10 5 Hz ( in Fig. 4 ). Comparing with recent 
experiment g/ir = 50MHz [31], the coupling is much 
more weak. Moreover, the ratio gj k is only 1.3, see Fig. 4, 
which can be realized in recent experiment technique [32j . 

III. CONCLUSION 

We proposed a scheme via three-level cascade atom to 
entangle two-mode fields as well as optomechanical os- 



cillator. Our study show that intracavity fields and two 
movable mirrors are entangled respectively for realizable 
coupling between the cavity fields and the atoms. The 
output two-mode fields entanglement is affected by the 
frequency of the mirror motion, the larger frequency of 
the mirror, the larger entanglement of the output fields. 
Because the entanglement is resulted from the atomic 
coherence, the initial atomic state play important role 
on the output entanglement of the two-mode fields. For 
the best coherence of the initial atomic state, the output 
spectrum appears one deep squeezing while for other ini- 
tial atomic state, the one deep squeezing split into two 
deep squeezing. 
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